Theory of stripe domains in magnetic shape memory alloys 
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The evolution of multivariant patterns in thin plates of magnetic shape memory materials with 
an applied magnetic field was studied theoretically. A geometrical domain-model is considered 
composed of straight stripe-like martensite variants with constant internal magnetization (high 
anisotropy limit) and magnetic domain wall orientation fixed by the twin boundaries. Through 
integral transforms of the demagnetization energy, the micromagnetic energy is cast into a form 
convenient for direct numerical evaluation and analytical calculations. The equilibrium geometrical 
parameters of multivariant patterns with straight and oblique twin boundaries have been derived as 
functions of the applied field and the material parameters of a plate. It is shown that the oblique 
multivariant states exist only in plates with thicknesses L larger than a certain critical value Lq. In 
samples with L < Lo a magnetic-field-driven transformation occurs directly between single variant 
states. 

PACS numbers: 75.60.Ch,75.60. Jk, 75.80.+q 



I. INTRODUCTION 

In ferromagnetic shape memory alloys like Ni-Mn-Ga 
Heusler alloys the high magnetocrystalline anisotropy 
fixes the magnetization along the easy magnetization axis 
of martensite variants [l| . 

This "one-to-one correspondence" between magnetic 
domains and the martensite variants allows to rearrange 
the crystallographic variants by applying a magnetic field 
[J, Hence, magnetic energy determines the evolu- 
tion of the multivariant states in magnetic fields. The 
martensitic microstructure can be described by adapt- 
ingthe phenomenological theory of magnetic domains 
[H,|i| . A detailed understanding of this coupling between 
microstructure and magnetic domain structure can be 
achieved for suitably simple geometrical systems. 

Here, we develop a theory of multivariant stripes in 
thin plates. For these stripe structure we derive analyt- 
ical expressions for demagnetization energies and com- 
plete phase diagrams for equilibrium structures. 



II. MICROMAGNETIC ENERGY AND 
EQUATIONS 

As a model we consider multivariant states in a layer 
of thickness L with surface normal along z and infinite 
extension in x and y direction. Commonly observed pat- 
terns include lamellar microstructures built from two co- 
existing single- variant states with oblique or straight in- 
terfaces, i.e. the angle between easy-magnetization axis 
and z is equal to or 7r/4, respectively (Fig. 1 a, b). 
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The orientation of these interfaces with 90-degree mag- 
netic walls are fixed by the crystallography of the twin- 
boundaries between the martensite variants. Within in- 
dividual variants stripe domains with 180-degree mag- 
netic domain walls can occur (Fig. 1, c) [1, IS S HI • 

The patterns in Fig. 1 a, b on the one hand sat- 
isfy the elastic compatibility between the tetragonal vari- 
ants and, on the other hand, comply with a common 
property of magnetic domains by avoiding uncompen- 
sated magnetostatic charges on the domain boundaries. 
Here, we assume that internal magnetic charges, that 
may arise at the twin-boundary Q , are absent or can be 
neglected as the period of 180-degree domain structures 
within variants remain small. Further, we assume that 
the uniaxial anisotropy is much stronger than the ap- 
plied fields. Thus, deviations of the magnetization from 
the easy axis within the variants can be neglected. As 
sketched in Fig. 1 a, the magnetization in the microsc- 
tructure with oblique twin interfaces can be reduced to 
a pattern with alternating magnetization m = ±(Af/2)z 
perpendicular to the plate surface and an effective bias 
field H = H z + 2irM . For the pattern with straight 
interfaces in Fig. 1 b, the alternating magnetization 
m = ±(M/V2)z, and the bias field is H = H z . Thus, 
the patterns in Fig. 1 are reduced to models of stripe do- 
mains with alternating magnetization m perpendicular 
to the surface and straight and oblique domain walls in a 
bias field H. The model with straight domain walls has 
been investigated in details for magnetic films (see e. g. 

and the solutions for the micromagnetic problems 
can be applied directly for magnetic shape memory films. 
The model with oblique interfaces is the main subject of 
our investigations in this paper. 

The reduced energy density w — WJ (2nm 2 ) for the 
oblique stripes can be written 
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FIG. 1: Multivariant patterns consisting of two different 
tetragonal martensites, (a) (3 = n/4, (b) (3 = 0, and stripe 
domains within one variant (c). 



where p = 2nL/D, q = (d + — tf_)/D, and D = d + + d- 
is the period of the domain structure, d± are the widths 
of domains polarized in the directions parallel (+) and 
antiparallcl (-) to the field. The domain wall energy (the 
first term in Eq. (fTJ) includes the characteristic length 
L c = <j/(Aitm 2 ) where a is the domain wall energy den- 
sity which describes the balance between the domain wall 
and stray field energies. The stray field energy Wd is given 
by a function S(p, g, 7r/4) that includes an infinite sum, 
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With the help of the integral transformation intro- 
duced in Q the infinite sum in Eq. is transformed 
into integrals on the interval [0,1]. Then, the energy ([I]) 
can be written 
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FIG. 2: Energy difference Aw(p) between stripe and single 
variant equilibrium states as function of the reduced inverse 
stripe period p and for different values of the reduced differ- 
ence q between widths of up and down domains. The equilib- 
rium energy of the multidomain states is lower than that of 
the homogeneous state. 
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1 — H I (47rm), can be written in the following form 
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where h = H/(Awm), L = 2\/2L Cl £ = pt, 

f(^l) — 2cos 2 (7rg/2) sinh£sin£x 

[ff(£> Q) ~ 2cos 2 (7r<7/2)(cosh£ cos£ + cos7rq)] _1 and 
q) = (cosh£ cos£ -I- cos7rg) 2 + sinh 2 £ sin 2 £. 

For straight stripes ( /3 = ) the integral transforma- 
tion of sum ([3]) the following expression gives the system 
energy [H 
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Minimization of Eqs. ((3]) or Q with respect to the inter- 
nal parameters p, q gives the equilibrium values for the 
domain sizes (d + , 



where G(t;,q) = (sinh £ — sin £ ) (cos 7rq + cosh £ cos f ) — 
sinh £ sin £ (sinh £ + sin £) . 

The energy difference Aw(p, q), Eq. ([5]), is negative for 
all p > and < < 1, and reaches zero only for p = 0, 
g = 1 (Fig. [2]). Similar relations are also true for straight 
stripes [9(. For both types of systems, stripes with finite 
widths have always lower energy than a single variant 
state. The stripes transform into a single variant state 
continuously by unlimited growth of the period D. Near 
the transition into the single variant state j? — > 0, and, 
thus, the function /(£) in Eq. ([3]) can be expanded with 
respect to the small parameter £, and the energy ([3]) can 
be evaluated in terms of elementary functions 
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III. RESULTS 

A. Transition into the single variant states 

The difference Aw = w — Wq between the system en- 
ergy w and that of the homogeneously magnetized layer, 
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u = 2 cos (nq/2)/p , and v(u) — 2(1 + it" 1 ) arctan(l + 
u^ 1 ) — 2(1 — u^ 1 ) arctan(l — u^ 1 ). Minimizing w from 
Eq. ^ with respect to u and q yields equations for equi- 
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librium values of p, q and u 
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Eqs. ([7]) express the equilibrium values of p and g in 
parametrized form as functions of the applied field h and 
the plate thickness L. According to Eqs. (|7J the param- 
eter h* from Eq. {2J defines the transition field into the 
saturated state (p = 0, q = 1). Note that for h = h* 
the parameter u transforms into u(h*) = d-/L. This 
means that the transtion into the saturated state takes 
place by an infinite expansion of the domains with the 
magnetization parallel to the field, (d + — > oo), while do- 
mains with the antiparallel magnetization (d—) keep a 
finite size. The transition line h* plotted as function of 
L c /L in Fig. [3] separates the regions of multivariant and 
single variant states. The equation h* = defines the 
critical thickness Lq = 2y2L c - For L < Lq the oblique 
domains do not exist, while straight domains (theoreti- 
cally) exist for any thickness Hi (Fig.©. 
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FIG. 3: Phase diagram in terms of the inverse reduced plate 
thickness L c /L and the reduced magnetic field H/(Anm) 
perpendicular to the plate for oblique (solid) and straight 
(dashed) stripes. The oblique domains do not exist below 
the critical thickness Lq = 2\/2L c . Inset shows the functions 
H*(L/L C ) close to the critical thickness. 
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FIG. 4: The equilibrium stripe domain sizes d+/L c (up- 
per set of curves) d—/L c (lower set of curves) as functions of 
the applied magnetic field for different layer thicknesses. In- 
set: dependence of the size for the minority domain on layer 
thickness at the transition. 



Typical solutions for domain sizes (d+ , cL ) are pre- 
sented in Fig. 2J At low fields, the domain sizes change 
only slowly. Exponential growth of the size d+ for the do- 
mains with magnetization in field direction sets in close 
to the saturation field. The minority domain size d- 
gradually decreases in the increasing magnetic field but 
remains finite at the transition field. The equilibrium 
magnetization curves in Fig. [5] display different behav- 
ior of the magnetization reversal in "thin" (L > L c ) and 
"thick" (L ^> L c ) plates. The former are characterized by 
low values of the saturation field /i*<Cl and steep mag- 
netization curves. In accordance with the phase theory 
the latter have linear magnetization curves m/m s = h in 
a broad range of the applied field. 

In various aspects, we find that the solutions for the 
oblique stripes have a universal character that may be 
useful for analysis of experimental data. The solutions 



for q in Eq. ([7]) give a universal magnetization curve for 
systems with large domains m/m s = (2/ir) arcsin(/i/ft,*) 
(Inset in Fig. [5]) . The solutions for the equilibrium values 
of p (8/V) Jp 1 1 arctan [sin ptj smhpt]dt = L /L at zero 
field ( q = 0) yield a specific dependence of zero-field 
period Dq on the layer thickness (Fig. [5]). 

The equation for p corresponding to the 
minimum value of Dq, arctan [sin pj sinhp] = 
Jp 1 1 arctan [sin ptj smhpt]dt, includes no material 
parameters and has the solution p m i n — 2.01. This 
yields L( mm )/L c = 4.58, D min) /L c = 14.33 and the 
ratio flj™"'^'"""' = 3.13. Corresponding values for 
straight domains are p min = 1.595, U min ^/L c = 3.8178, 
0™ in) /L c = 14.5739, = 3.8176 

[lfj |. The minority domain size , d*_/L c (L/L c ), 
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FIG. 5: Magnetization curves for different values of the re- 
duced thickness L/L c . Inset shows the "universal" character 
of the magnetization curves for large domains. 



at the transition field has a similar minimum (In- 
set in Fig. 2]). The equation for x — vam{d*_/L) 
4arctan (x 2 /2) — x 2 In (l + 4/x 4 ) = has the solution 
x = 0.918, and yields the following values of the 
parameters in this point: L/L c = 5.573, d*_/L c = 5.115. 

The point marks the transition region between 

a functional dependence Dq oc sqrtL characteristic for 
thick plates (Do -C L) (Kittel law) to a exponential 
growth of the period with a decreasing plate thickness, 
Do oc exp(—7ri c /L)Q in thin layers (L > L c ). This tran- 
sitional region has been reached in thin polycrystallinc 
films of Ni51.4Mn28.3Ga— 20.3 with easy-axes dominantly 
at 45-degree from the layer normal direction [3] ■ Inset in 
Fig. [5] shows the calculated functions Dq(L) with mate- 
rial parameter L c adjusted to the experimental data of 
Ref. Q. 



IV. CONCLUSIONS 




The present theory enables rigorous calculations of 
equilibrium structures in thin ferromagnetic martensite 
plates with two variants. It is clear that real magnetiza- 
tion processes in magnetic shape-memory processes will 
generally be hysteretic. However, the equilibrium states 
could be reached, e.g., free-standing films or thin free 
single-crystal platelets with very mobile twin-boundaries. 
The analytical expressions for the energy of a stripe sys- 
tem Eqs. ((T])-(energyST) are applicable also for non- 
equilibrium states. We find that oblique stripe structures 
are possible only above a certain critical thickness, while 
straight stripe structures have been known to exist for ar- 
bitrary dimensions of layers with perpendicular magnetic 
anisotropy. The theoretical method can be extended to 
treat twinned martensite plates with other angle between 
easy magnetization and layer normal direction. 



FIG. 6: The equilibrium period Do/L c as a function of the 
plate thickness L/L c at zero field (solid line) in comparison 
with a correponding line (dashed) for straight stripes. Inset 
shows D(L) functions for a thin layer of Nisi. 4Mn28.3Ga— 20.3 
investigated in Ref. |7| corresponding to the straight stripe 
structure Fig. [Tj b. The dashed curve is plotted for the av- 
eraged characteristic length L c , the solid line incorporates 
effects connected with the indicated variation of L c with L. 



Acknowledgements 

The authors are grateful to J. McCord, Y. W. Lai, 
V. Neu, and R. Schafcr for helpful discussions. N.S.K., 
A.N.B. thank H. Eschrig for hospitality at IFW Dres- 
den. We gratefully acknowledge support by Deutsche 
Forschungsgmeinschaft (SPP 1239 project A8). 



[1] Q. Pan, R. D. James, J. Appl. Phys. 87, 4702 (2000). 
[2] O. Heczko, K. Jurek, K.Ullakko, J.Magn. Magn. Mater. 

226, 996 (2001). 
[3] A. N. Bogdanov, A. DeSimone, S. Miiller, U. K. Rofiler, 

J. Magn. Magn. Matter. 261, 204 (2003). 
[4] A. Hubert, R. Schafer, Magnetic Domains (Springer- 

Verlag, Berlin, 1998). 
[5] Y. Murakami et al., Appl. Phys. Lett. 82, 3695 (2003). 
[6] Y.Ge, O. Heczko, O.Soderberg, and V.K. Lindroos, 



J.Appl.Phys. 96, 2159 (2004). 
[7] V.A. Chernenko et al., Acta Materialia 54, 5461 (2006). 
[8] Y.W. Lai et al., Appl. Phys. Lett. 90, 192504 (2007). 
[9] A. N. Bogdanov, D. A. Yablonskii, Fiz. Tverd. Tela 22, 

680 (1980) [Sov. Phys. Solid State 22, 399 (1980)]. 
[10] A. N. Bogdanov, U. K. Rofiler, |cond-mat/0 606671 

(2006). 

[11] N. S. Kiselev, I. E. Dragunov, V. Neu, U. K. Rofiler, A. 
N. Bogdanov, to be published. 



